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1.  Statement  of  the  problem  studied 

Studies  of  sound  propagation  in  a  turbulent  atmosphere  are  important  for  both  military  and  civil 
applications.  First,  these  studies  can  significantly  improve  performance  of  U.S.  Army  acoustic  sensors 
for  source  detection,  ranging  and  recognition.  Furthermore,  many  U.S.  Army  ’’smart”  weapons  rely  on 
acoustic  sensors.  Second,  civil  applications  of  studies  of  sound  propagation  in  a  turbulent  atmosphere 
include  noise  pollution  near  highways  and  airports,  acoustic  remote  sensing  of  the  atmosphere,  and 
theories  of  waves  in  random  media. 

Most  of  the  previous  research  in  atmospheric  acoustics  has  been  devoted  to  studies  of  sound 
propagation  in  daytime  boundary  layers.  In  our  two  previous  ARO  projects,  we  developed  theories 
of  sound  propagation  in  daytime  boundary  layers  [see,  G.H.  Goedecke  and  V.E.  Ostashev,  ARO 
Project  D A AH04-95- 1-0593  (1995-1998);  V.E.  Ostashev  and  G.H.  Goedecke,  ARO  Project  DAAG55- 
98-1-0463  (1998-2001)].  However,  studies  of  sound  propagation  and  scattering  in  nighttime  boundary 
layers  (NBL)  are  also  very  important  because  Army  sensors  and  smart  weapons  must  operate  during 
day  and  night.  Note  that  nighttime  conditions  are  usually  very  favorable  to  acoustic  detection.  The 
main  gaols  of  this  grant  were  development  of  3D  models  of  turbulence  and  internal  gravity  waves 
(IGW)  in  NBL  and  studies  of  sound  propagation  and  scattering  in  NBL. 

Specific  tasks  of  the  grant  and  their  accomplishment  are  presented  below  in  section  2. 


2.  Summary  of  the  most  important  results 

As  a  result  of  accomplishment  of  this  grant,  we  developed  a  theory  of  sound  propagation  and  scattering 
in  NBL. 

First,  we  obtained  3D  spectra  of  temperature  and  wind  velocity  fluctuations  due  to  IGW  in  mod¬ 
erately  and  very  stable  NBL.  Also,  we  developed  a  3D  model  of  isotropic  inhomogeneous  turbulence  in 
weakly  stable  NBL.  Finally,  we  developed  quasi-wavelet  (QW)  models  of  turbulent  temperature  and 
wind  velocity  fluctuations. 

Then,  the  developed  models  of  IGW,  turbulence,  and  quasi-wavelets  (QWs)  were  used  to  study 
sound  propagation  and  scattering.  We  considered  three  different  problems:  scattering  of  sound,  line- 
of-sight  sound  propagation,  and  sound  propagation  in  a  refractive,  turbulent  atmosphere  near  an 
impedance  ground.  Furthermore,  the  coherence  function  of  a  spherical  sound  wave  propagating  in 
NBL  was  measured  and  compared  with  that  predicted  theoretically. 

The  specific  tasks  of  the  grant  and  their  accomplishment  are  as  follows: 

Task  1.  3D  models  of  temperature  and  wind  velocity  fluctuations  due  to  IGW  in  moderately  and 
very  stable  NBL. 

Accomplishment.  This  task  was  accomplished  in  Refs.  [24-27].  First,  we  obtained  the  3D  spectra 
of  temperature  and  wind  velocity  fluctuations  due  to  IGW  [24-26].  The  normalized  3D  spectrum  of 
temperature  fluctuations  $t(A_l,  A3)  is  shown  in  Fig.  1.  Here,  K%  and  K±  are  the  components  of  the 
turbulence  wave  vector  in  the  vertical  and  horizontal  directions,  respectively.  In  Fig.  1,  Kq  is  inversely 
proportional  to  the  outer  length-scale  of  the  spectrum.  It  follows  from  the  figure  that  K%) 

reaches  the  maximum  at  K±  =  K%  =  0  and  decreases  along  the  A^-axis  much  faster  than  along 
the  A3-axis;  i.e.  the  spectrum  is  highly  anisotropic.  Second,  we  developed  a  model  which  predicts 
fluctuations  in  temperature  and  wind  velocity  due  to  one  or  several  IGW  in  stably  stratified  NBL  [27]. 
The  model  can  be  used  for  calculations  of  sound  propagation  through  a  snapshot  of  temperature  and 
wind  velocity  fluctuations  due  to  IGW. 

Task  2.  3D  models  of  isotropic,  inhomogeneous,  and  intermittent  turbulence  in  weakly  stable  NBL. 
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Figure  1:  3D  spectrum  of  temperature  fluctuations  due  to  IGW. 

Accomplishment.  In  weakly  stable  NBL,  the  3D  spectra  of  temperature  and  wind  velocity  fluc¬ 
tuations  were  modeled  as  isotropic,  inhomogeneous  von  Karman  spectra,  see  Eqs.  (l)-(8)  in  [3]  or 
Eqs.  (l)-(4)  in  [18].  In  these  equations,  for  the  considered  case  of  a  weakly  stable  stratification,  one 
should  ignore  terms  containing  the  Monin-Obukhov  scale  La  in  the  temperature  spectrum  and  retain 
only  shear-produced  velocity  fluctuations.  The  use  of  the  von  Karman  spectra  of  temperature  and 
wind  velocity  fluctuations  allows  us  to  correctly  describe  the  inertial  subrange  of  turbulence  and  to 
realistically  account  for  inhomogeneities  in  the  energy  subrange.  Figure  2  from  Ref.  [3]  depicts  the 
von  Karman  spectrum  of  temperature  fluctuations  (referred  to  as  ’’simplified  von  Karman”)  and  a 
more  general  von  Karman  spectrum  of  temperature  fluctuations  which  has  a  Gaussian  cutoff  at  large 
wavenumbers.  The  von  Karman  spectrum  of  velocity  fluctuations  is  shown  in  Fig.  3  obtained  in 
Ref.  [23]. 

It  is  well  known  that  the  turbulence  in  NBL  can  be  intermittent,  i.e.  it  can  occur  in  bursts  of 
activity.  In  [14,17],  the  intermittency  of  turbulence  was  modelled  by  assuming  that  the  structure 
parameters  of  temperature  and  wind  velocity  fluctuations,  Cj.  and  C%,  are  random  functions.  Note 
that  the  considered  above  von  Karman  spectra  of  temperature  and  wind  velocity  fluctuations  are 
proportional  to  Cf-  and  C%,  respectively.  Therefore,  these  spectra  with  random  values  of  Cf-  and  C% 
provide  3D  spectra  of  temperature  and  wind  velocity  fluctuations  in  weakly  stable,  intermittent  NBL. 

Task  3.  Quasi-wavelet  models  of  turbulent  temperature  and  wind  velocity  fluctuations. 
Accomplishment.  The  task  was  completed  in  Refs.  [3,7,16,18,22,23].  During  the  reporting  period,  we 
have  continued  the  development  of  QW  models  of  turbulence  which  was  started  in  our  previous  ARO 
project.  QWs  are  localized,  self-similar  temperature  or  wind  velocity  eddies.  The  ensemble  of  QWs 
with  many  different  sizes  and  properly  chosen  scaling  laws  resembles  actual  turbulence  [3,7,16,18,23]. 
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normalized  temperature  spectra 


Figure  2:  Spectra  of  temperature  fluctuations. 


Figure  3:  Spectra  of  velocity  fluctuations. 
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The  use  of  QW  models  in  numerical  and  analytical  calculations  of  sound  propagation  and  scattering 
in  a  turbulent  atmosphere  can  significantly  simplify  such  calculations  and  be  advantageous  in  compar¬ 
ison  with  other  methods.  During  the  project,  we  showed  that  isotropic  turbulence  with  any  physically 
’’reasonable”  spectra  of  temperature  and  wind  velocity  fluctuations  can  be  modeled  as  an  ensemble 
of  QWs  with  specially  chosen  parent  functions  which  determine  the  shape  of  QWs.  In  particular,  we 
found  the  parent  functions  and,  hence,  QW  models,  which  exactly  result  in  the  von  Karman  spectra  of 
temperature  and  wind  velocity  fluctuations,  see  Refs.  [3,16,22,23].  This  is  an  important  result  which 
can  be  considered  as  a  phenomenological  motivation  for  the  von  Karman  spectra.  Furthermore,  we 
found  simple  parent  functions,  gaussian  and  exponential  (in  the  Fourier  space),  which  yield  the  spectra 
that  are  very  close  to  the  von  Karman  spectra,  see  Figs.  2  and  3.  Finally,  in  [7],  nonspherical  QWs 
were  used  to  model  anisotropic,  shear-driven  velocity  fluctuations.  Note  that  wind  shear  is  one  of  the 
mechanisms  of  turbulence  in  NBL. 

Task  4.  Sound  scattering  in  a  turbulent  atmosphere. 

Accomplishment.  The  task  was  completed  in  Refs.  [1,2,8,10,12,19-21,24-26].  First,  we  derived  the 
exact  formula  for  the  sound  scattering  cross-section  per  unit  volume,  cr,  in  a  turbulent  atmosphere 
[8,20] .  This  formula  was  derived  from  complete  sets  of  linearized  equations  of  fluid  dynamics  considered 
in  [1,12,21].  In  comparison  with  equations  for  a  known  in  the  literature,  the  derived  formula  contains 
several  new  terms  which  describe  sound  scattering  by  atmospheric  pressure  fluctuations,  the  potential 
component  of  velocity,  the  cross-correlation  between  temperature  and  wind  velocity  fluctuations,  etc. 
Using  this  formula,  we  studied  sound  scattering  by  atmospheric  pressure  fluctuations  in  NBL  [8]. 

Second,  we  calculated  the  sound  scattering  cross-section  a  for  the  case  of  sound  scattering  by 
temperature  and  wind  velocity  fluctuations  due  to  IGW  [24-26].  Figure  4  depicts  the  normalized  a  as 
a  function  of  9  —  6q.  Here,  9q  =  85°  and  6  are  the  angles  between  the  vertical  axis  and  the  directions 
of  propagation  of  the  incident  and  scattered  waves,  respectively.  A  sharp  maximum  in  a  occurring 
at  9  —  9q  =  0  is  a  well  known  result  in  theories  of  waves  in  random  media.  In  Fig.  4,  there  is  also  a 
second,  much  broader  maximum  which  occurs  at  9  —  9o  =  10°.  This  maximum  is  a  new,  interesting 
result  which  does  not  appear  for  the  case  of  sound  scattering  by  isotropic  turbulence. 

Third,  the  QW  model  of  turbulent  velocity  fluctuations  was  used  to  study  sound  scattering  behind 
barriers  [2,19],  and  sound  scattering  into  refractive  shadow  zone  [10]. 

Task  5.  Line-of-sight  sound  propagation. 

Accomplishment.  References  [4,5,14,15,17,24-26]  accomplish  this  task.  For  line-of-sight  sound  prop¬ 
agation,  the  most  widely  used  statistical  moments  of  a  sound  field  are  the  mean  sound  field,  the  coher¬ 
ence  function,  and  the  structure  functions  of  log-amplitude  and  phase  fluctuations.  In  Refs.  [4,5,14], 
general  formulas  for  these  statistical  moments  of  plane  and  spherical  sound  waves  propagating  through 
inhomogeneous,  anisotropic  turbulence  with  arbitrary  spectra  of  temperature  and  wind  velocity  fluctu¬ 
ations  were  derived.  The  derived  general  formulas  were  then  used  to  calculate  and  study  the  statistical 
moments  of  a  sound  field  for  some  important  particular  cases. 

First,  using  these  formulas,  the  mean  sound  field  and  the  coherence  function  were  calculated  for 
the  inhomogeneous  von  Karman  spectra  of  temperature  and  wind  velocity  fluctuations  [4]  and  Mann’s 
anisotropic  spectral  tensor  of  shear  driven  turbulence  [14], 

Second,  the  general  formulas  for  the  above-mentioned  statistical  moments  of  plane  and  spherical 
sound  waves  were  used  to  calculate  these  statistical  moments  for  the  case  of  sound  propagation  through 
intermittent  turbulence  [14,17].  The  calculations  were  done  by  assuming  that  the  structure  function 
parameters  and  C%  are  random  functions  along  the  sound  propagation  path  as  described  in  Task  2. 

Third,  the  general  formulas  were  used  to  calculate  the  mean  sound  field  and  the  coherence  function 
of  a  plane  sound  wave  propagating  in  NBL  with  IGW  [24-26].  Figure  5  from  Ref.  [24]  shows  the 
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Figure  4:  Normalized  scattering  cross-section  due  to  sound  scattering  by  IGW. 


normalized  coherence  function  T(x\y,z)  versus  the  normalized  separation  between  sensors.  In  figure, 
r(x;  y,  0)  is  the  normalized  coherence  function  for  the  horizontal  separation  between  two  sensors  along 
the  y-axis,  and  r(x;  0,  z)  is  that  for  the  vertical  separation  between  sensors  along  the  z-axis.  Different 
values  of  eo  correspond  to  a  different  degree  of  anisotropy  of  IGW  field.  It  follows  from  Fig.  5  that  the 
coherence  decays  much  more  rapidly  with  increasing  vertical  separation  between  sensors  than  with 
increasing  horizontal  separation. 

Finally,  note  that  the  state-of-the-art  array  for  acoustic  tomography  of  the  atmosphere  which 
is  under  construction  by  several  organizations  in  the  U.S.  (ARO  projects  DAAD19-03-1-0104  and 
DAAD19-03-1-0341)  will  be  based  on  measurements  of  travel  times  of  line-of-sight  sound  propagation 
between  different  pairs  of  sources  and  receivers.  Different  theoretical  aspects  of  this  array  were  studied 
in  [15]. 

Task  6.  Combined  effects  of  refraction,  turbulence,  and  impedance  ground  on  sound  propagation. 
Accomplishment.  This  task  was  accomplished  in  Refs.  [1,6,9,11-13,21].  First,  using  some  me¬ 
teorological  parameters,  classification  schemes  for  sound  propagation  in  the  nighttime  and  daytime 
atmosphere  near  the  ground  have  been  developed  [6,11].  The  transmission  loss  of  sound  propagation 
was  classified  in  terms  of  these  parameters. 

Second,  fast,  accurate  predictions  of  the  sound  pressure  level  in  a  refractive,  turbulent  atmosphere 
were  developed  based  on  an  artificial  neural  network  [9].  The  neural  network  training  parameters 
were  the  following:  source  and  receiver  height,  the  horizontal  distance  between  source  and  receiver, 
the  azimuthal  direction  of  sound  propagation,  the  static  flow  resistivity  of  the  ground,  the  friction 
velocity,  the  turbulent  sound  speed  scale,  and  a  range  of  frequencies. 

Third,  new  starting  equations  for  calculations  of  sound  propagation  in  a  refractive,  turbulent, 
and  moving  atmosphere  were  derived:  a  wide  angle  parabolic  equation  [13]  and  equations  for  finite- 
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Figure  5:  Normalized  coherence  function  of  a  plane  sound  wave  propagating  through  IGW  field. 

difference  time  domain  calculations  (FDTD)  [1,12,21].  Note  that  FDTD  techniques,  which  can  provide 
the  most  accurate  calculations  of  the  sound  field  in  complicated  environments,  have  been  under  in¬ 
tensive  development  by  Army  scientists. 

Task  7.  Comparison  of  experimental  results  with  theoretical  predictions. 

Accomplishment.  This  task  was  completed  in  Ref.  [28].  First,  the  coherence  function  of  a  spherical 
sound  wave  and  the  spectra  of  temperature  and  wind  velocity  fluctuations  were  measured  during  24 
hours  or  even  longer.  The  experimental  results  obtained  were  compared  with  theoretical  predictions 
of  the  coherence  function  obtained  in  task  5. 


Technology  transfer 

In  accomplishment  of  this  grant,  we  closely  collaborated  with  the  U.S.  Army  scientists:  Dr.  K.  Wilson 
(CRREL),  Dr.  M.  Mungiole  (ARO),  D.  Marlin  (ARO),  and  H.  Auvermann  (ARO).  This  means  a  direct 
transfer  of  the  results  obtained  for  needs  of  the  CRREL  and  ARL.  The  Army  scientists  are  coauthors  in 
25  of  the  total  of  28  papers  devoted  to  accomplishment  of  the  grant,  see  the  List  of  Publications  below. 
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